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Abstract. — In this work, we prove global existence of solutions for second order differential 
problems in a general framework. More precisely, we consider second order differential inclusions 
involving proximal normal cone to a set- valued map. This set- valued map is supposed to take 
admissible values (so in particular uniformly prox-regular values, which may be non-smooth and 
non-convex). Moreover we require the solution to satisfy an impact law, appearing in the description 
of mechanical systems with inelastic shocks. 

Resume. — Nous prouvons dans ce travail un resultat d'existence globale de solutions pour 
des inclusions differentielles du second ordre. Nous considerons plus precisement des inclusions 
differentielles du second ordre faisant apparaitre le cone proximal normal d'un ensemble C 
dependant du temps. La multifonction C(-) est supposee admissible (et en particulier prendre des 
valeurs uniformement prox-regulieres eventuellement non regulieres et non convexes). De plus, nous 
imposons a la solution de verifier une lot d'impact, apparaissant dans la description de systemes 
mecaniques avec des chocs inelastiques. 
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1. Introduction 

We consider second order differential inclusions, involving proximal normal cones. These 
differential inclusions appear in several fields (granular media [25, 27], robotics [16], [17] and 
virtual reality [34] ... ). 

They describe an evolution problem where the state- variable is submitted to some constraints 
and so has to live in an admissible set). 
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Numerous works deal with the particular case where the admissible set is given by several 
constraints. Let us first detail them. We write C C M. d for the closed admissible set defined by 
constraints ( 9 i)i<i<p as follows : 



v 

C : = p| {q G R d , 9i (q) > o} . (1) 



=1 

The tangent cone to C at q is 

T c (q) := {u G M d , (u, V 9i {q)) >0 Vie /(<?)} 

where I(q) is the set of "active constraints" 

/(g) := {i, 9i (q) = 0}. 

The second order problem is the following one: let I be a bounded time-interval, / ilxl" 1 -} M. d 
be a map, find q G W 1 ' 00 ^, K d ) such that q G and 

( dq + N c (q)dt B f(t,q)dt 

VtGX, q(t + ) = P Tc{q{t)) (q(t")) 

q(0) = q G Int(C) 

k g(0) = u , 

where Int(C) is the interior of the set C and Nc(q) is the normal cone defined by 

( {0} if q G Int(C) 

N c (q) := { {- £i 6 J(«) A i V #i(9), A * > °} if 9 € 3C 

if g £ C. 

This differential inclusion can be thought as follows: the point q(t), submitted to the external 
force f(t, q(t)), has to live in the set C. The differential inclusion dq + Nc{q)dt B f(t, q)dt does 
not uniquely define the evolution of the velocity during an impact. To complete the description, 
we impose the impact law q(t + ) = PT c (q(f))(?(t _ ), introduced by J.J. Moreau in [25] and justified 
by L. Paoli and M. Schatzman in [28, 30] (using a penalty method) for inelastic impacts. This 
law can be extended with a restitution coefficient to model the elastic shocks. 

The existence of a solution for such second-order problems is still open in a general framework. 
The first positive results were obtained by M.P.D. Monteiro Marques [22], L. Paoli and M. 
Schatzman [29] in the case of a smooth admissible set (which locally corresponds to the single 
constraint case p = 1 in (1)). The single constraint case is also treated in [23, 13] where an 
additional mass-matrix depending on q. The proofs use a time-discretization of (2) and rely on 
the convergence of the approximate solutions. The multi-constraint case with analytical data 
was then treated by P. Ballard with a different method in [2]. In this paper, a positive result 
of uniqueness for such problems was also obtained. Then in [31], an existence result is proved 
in the case of a non-smooth convex admissible set C (given by multiple constraints). There, 
the active constraints are supposed to be linearly independent in the following sense: for each 
q G dC, the gradients (V 'gi(q))i^i( q ) are supposed to be linearly independent. This assumption 
is quite strong since it implies that the number of active constraints \I(q)\ is always lower than 
the dimension d (which may fail for some applications). Moreover the impact law is proved 
under a geometrical assumption: 

V(t, j) G I(q)\ i + j, (V 9i (q), V 9j (q)) < 0. (3) 

In [32, 33], similar results are obtained without requiring the convexity of C. Note that excepted 
in [2], all these results are local in the following sense: for each initial data (qo,uo), there exists 
a time r = r(|ito|) and a solution q of (2) on I = [0, r]. 
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Recently, time-dependent constraints are considered. More precisely, the constraints gi 
X x M. d — > R define the following set-valued map 



v 

C(-):=f]{q€R d , 9i(;q)>0}. (4) 
i=l 

Then Problem (2) becomes: find q G W^°°{l,R d ) such that q G BV(l,R d ) 

' dq + N c[t) {q)dt B f(t,q)dt 
WgZ, q (t+)=P VtqW q(t-) 
9(0) = go e Int(C(0)) 
9(0) = u Q , 

where Vt t q is the set of admissible velocities: 

V t>q := {u, d t9i {t, q) + (V, gi (t, q),u}>0 if q) = 0} . (6) 

Local existence results were obtained in [40] assuming that dC(-) belongs to C 3 (l x R d ) (which 
locally corresponds to a single constraint). Recently in [5], a global result was proved by the 
first author and A. Lefebvre for convex C 2 constraints <?, (no regularity on the boundaries dC{t) 
is required). Moreover, Assumption (3) and the independence of the gradients (y q gi(t, q))iei(t,q) 
are relaxed to a positive linearly independence: global results are shown under the existence of 
p, 7 > such that for all q G C(t) and nonnegative reals Aj 



^2 XiV q gi(t,q) 

iel p (t,q) 



MV q 9i(t,q)\<i 

iel p (t,q) 

where 

I p (t,q) := 0, 9 i(t,q) < p} . 



In this paper, we are interested in local and global existence results for second order problems 
involving a general set-valued map C(-) (not necessary defined with constraints). This problem 
was solved in the case of first order differential inclusions called sweeping process. Let us briefly 
recall the corresponding context. For a bounded time- interval X, a set- valued map C : X R d 
with nonempty closed values and a map / ilxl^ R d , the associated sweeping process takes 
the following form: 

§(*)+N(C(t),g(t)) 9 /(*,(?(*)) 
9(0) = <?o , 

with an initial data 90 G C(0) and where N(C, 9) denotes the proximal normal cone of C at any 
point 9. This kind of evolution problem has been introduced by J.J. Moreau in 70's (see [26]) 
with convex sets C(t). He proposed a time-discretization of (7) called the catching-up algorithm. 
This scheme was later adapted to the second order problems (2) and (5). Later the convexity 
assumption was weakened by the concept of "uniform prox-regularity" (a set C is said to be 
uniformly prox-regular with constant rj or rj-prox-regular if the projection onto C is single- valued 
and continuous at any point distant at most rj from C). In this framework, the well-posedness of 
(7) was proved by G. Colombo, V.V. Goncharov in [11], H. Benabdellah in [3], L. Thibault [41] 
and by G. Colombo, M.D. P. Monteiro Marques in [12]. The sweeping process problem is still 
extensively studied but the recent results are not detailed here (see the works of M. Bounkhel, 
J.F. Edmond and L. Thibault in [7, 41, 14, 15] and of the authors in [4, 6]). 

Note that a major difference exists between the first and the second order differential inclusions. 
Indeed, even with smooth data the uniqueness does not hold for second order problem such (2) 
or (5), see [38] and [2]. The only positive results are proved in [39] for one-dimensional impact 
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problems and in [2] in the context of analytic data. Here we are only interested in existence 
results for general second order problems. 

In the case of a general set- valued map C(-), we introduce the following set (which extends the 
concept of Tc(q) and Vt tq ): 

o 

Definition 1.1. — For a set-valued map C : X z4 M. d , for every to £ X and go £ C(i), t/ie cone 
of admissible velocities is defined as follows: 

P J r +h r O (tp + e)-q \ . C(t + e)-q Q 

tt , qo := <v = hmv e , with v e £ > = iimmt . 

[ e\o e J e\o e 

We assumed that the set-valued map C(-) is Lipschitz continuous on X: there exists Co > such 
that for all t, s £ X 

d H {C(t),C(a)) <c |i-s|, (8) 
where d# denotes the Hausdorff distance. 

We are interested in the following problem: find q £ W 1,oc (I, M. d ) such that q £ BV{X, W 1 ) such 
that 

( dq + N(C(t),q)3 f(t,q)dt 

Vt£X, <K*+) = Pe-^r)] (Q) 
?(0) = go, 
k q(0) = u 

where go £ Int[C(0)] and uo are initial data and N(C(t), g) is the proximal normal cone to C(t) 
at q (see Definition 2.3). Let us give a more precise sense to this differential problem. 

Definition 1.2. — LetX := [0,T] 6e a bounded time-interval. A continuous function q : X — >• 
M d is a solution of (9) if there exists another function k : X — > M. d such that : 

a) q belongs to W 1,00 (X,M ,i ) 

b) q and k belong to BV(l,R d ) 

c) the following differential equation is satisfied in the sense of time-measures 

dq + dk = f(t,q)dt (10) 

o 

d) for all t £ X, the impact law q(t + ) = q ^[q(t~)] holds 

e) the initial conditions are verified : q(0) = qo and q(0) = uq 

f) the differential measure dk is supported on {t, q(t) £ dC(t)}: 

\k\(t) = I l gis)e ac(s)d\k\(s), k(t) = [ £(s)d\k\(s), (11) 
Jo Jo 

where £ : X — > R d is a measurable function satisfying for all s £ X: £ N(C(s), g(s)), 
|£(s)| = 1 and |fc|(i) := 7ar(fc, [0,t]). 

The main subject of this work is to prove the following global existence result. 

Theorem 1.3. — Let X = [0,T] be a bounded time interval. Suppose that C : X M. d is a 
Lipschitz and admissible (see Definition 2.13) set-valued map Assume that f : X xl d ->l ,i is 
a measurable map satisfying: 

3K L >0,Vtel,Vq,qeC(t), \f(t,q)-f(t,q)\<K L \q-q\ (12) 

3F£L 1 (X), W£X, Vg€C(t), \f(t,q)\ < F(t). (13) 

Then the differential inclusion (9) admits at least one solution. 
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The paper is structured as follows: the notions of uniform prox-regularity and admissibility 
related to the main assumptions of Theorem 1.3 are detailed in Section 2. Moreover, we prove a 
natural expression of the set Ct 0t q with the help of the tangent cone of Q, := {(t, q) S R x R d , t S 
I, q £ C(t)} at (to, <Zo)- Then Section 3 is devoted to the proof of Theorem 1.3. In Section 4, we 
deal with a set-valued map C(-), defined by constraints as in (4). More precisely, we prove under 
(Rp) global existence results (Theorem 4.6). For time-independent constraints, Assumption (Ro) 
is shown to be sufficient to obtain local existence results (see Theorem 4.7), which is weaker 
than the usual geometrical assumptions. 



2. Preliminaries and Definitions 

First we precise some notations. For a time interval I, we write W 1,oc (X, R d ) (resp. W l,1 (I, R )) 
for the Sobolev space of functions in L°°(Z, R d ) (resp. L l (T,M. d )) whose derivative is also in 
L°°(l,R d ) (resp. L l (X, R d )). BV{Z,R d ) is the space of functions in L°°(l,R d ) with bounded 
variations on I. We define the dual space M(I) '■= (C C (I))' where C C (I) is the space of contin- 
uous functions with compact support (corresponding to the set of Radon measure due to Riesz 
Theorem). We set Ai+(Z) for the subset of positive measures. 

We emphasize that the different notions defined in this section can be extended in the case of 
an infinite dimensional Hilbert space H. We consider the Euclidean space R d , equipped with 
its euclidean metric | | and its inner product {■,■}■ For a subset Q of R d , we denote by cIq the 
distance function to this set : 

d Q (x) := inf \y - x\. 

y&Q 

Definition 2.1. — Given a family of sets D e C R d indexed by e > 0, the outer and inner limits 
are defined respectively by 

limsupDe '■= \ x S R d , 3xk — > x, 3efc \ with Xk S D ek > , 

and 

liminf D f := \ x = lim x f , x f € D f > . 

If these two sets both equal a set D C R d , then we say that D e converges to D and we write 
lim D e = D. 

e\0 

Definition 2.2. — Let Q be a closed subset ofR d . The set-valued projection operator Pq is 
defined on R d by 

Vx e R d , P Q (x) := {y G Q, \x - y\ = d Q (x)} . 

Definition 2.3. — Let Q be a closed subset ofR d and x € Q, we write N(Q, x) for the proximal 
normal cone of Q at x, defined by: 

N(Q,x) := j?;e]R a! , 3s > 0, x e Pq(x + sv)}. 

The cone N(Q, x) somehow generalizes the notion of the outward normal direction. In Figure 1, 
we have plotted a set Q with several points X{. At the regular point X4 (where the boundary is 
smooth), the proximal normal cone is exactly the half-line directed by the outward direction. At 
the points xo, x\ and X3 the boundary is not smooth and the proximal normal vectors constitute 
a cone. At the point X2, notice that the proximal normal cone is reduced to {0}. 



We now define the tangent cone as follows: 
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N(Q,Xi) 




Figure 1 . Examples of proximal normal cones 



Definition 2-4- — Let Q be a closed subset ofM d and x G Q, we write Tq(x) for the tangential 
cone of Q at x, defined by the following outer limit: 

Q — x ( Q — x] 

Tq(x) : = limsup = < v G R d , 3v k ->■ v, 3e k \ 0, v k £ > . 

e\0 el e k J 



2.1. Uniform Prox-regularity. — We now come to the main notion of "prox-regularity" . It 
was initially introduced by H. Federer (in [18]) in finite dimensional spaces under the name of 
"positively reached sets". Then it was extended in infinite dimensional space and studied by 
F.H. Clarke, R.J. Stern and P.R. Wolenski in [8] and by R.A. Poliquin, R.T. Rockafellar and 
L. Thibault in [35]. 

Definition 2.5. — Let Q be a closed subset o/M d and rj > 0. The set Q is said rj-prox-regular 
if for all x G Q and v G N(Q,x) \ {0} 

B (x + rjj^,'i]\ (~)Q = 0. 

By extension, a closed set Q is said uniformly prox-regular if there exists r] > such that Q is 
■q-prox-regular. 



We refer the reader to [8, 9] for other equivalent definitions related to the limiting normal cone. 
The previous definition is very geometric, it describes the fact that we can continuously roll an 
external ball of radius rj on the whole boundary of the set Q. The main property is the following 
one: for a ry-prox-regular set Q, and every x satisfying (Iq(x) < rj, the projection of x onto Q is 
well-defined (i.e. Pq(x) is a singleton) and the projection is continuous. 

The notion of prox-regularity can be described by the hypomonotonicity property of the proximal 
normal cone too (see the work of R.A. Poliquin, R.T. Rockafellar and L. Thibault ([35])): 

Proposition 2.6 ([35]). — Let Q be a closed subset ofM. d . Then Q is rj -prox-regular if and 
only if the proximal normal cone N(Q, •) is an hypomonotone operator with constant j-. That 
means: for all x,y G Q and all (ce,/3) G N(Q, x) x N(Q,y), we have 

(a - 0,x - y) > -— [\a\ + \/3\] \x - y\ 2 
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or equivalently 



(a,V - x) < —\a\\x - y\ 2 . (14) 



Proposition 2.7. — Let Q be a n -prox-regular set ofM. d . Consider y £ W 1 with < 2n. 

Then y G Q if and only if 

\/x G QnB(y,d Q (y)), Va G N(Q,x), (a,y - x) < — \a\\x - y\ 2 . 

Proof: Consider a point y ^ Q satisfying the above property and choose x G Pg(y). Then 
y — x G N(Q, x) and we get 

i 1 2 ^ i i ^ 

\y — x\ < — \x — v \ , 
iv i - 2rj \ y\ i 

which leads to a contradiction with the assumption dq(y) < 2rj. □ 

Since the prox-regularity of a set implies its Clarke's regularity (see Definition 6.4 of [36]), we 
have the following lemma (Corollary 6.29 of [36]). 

Lemma 2.8. — Let Q be a uniformly prox-regular set. Then for all x G Q the normal and 
tangential cones are mutually polar : 

N(Q,x)°=T Q (x). 

That means: 

T Q (x) = |« € R d , W G N(Q,a?) (u,v) < o} . 

In particular, it comes that every tangential cone of a uniformly prox-regular set is convex. 

Proposition 2.9. — Let C(-) be a Lipschitz continuous set-valued map taking uniformly prox- 
regular values on 1. For all to £ I and qo G C(to), we set 

r . . C(t + h)-q - C (t + h)-q 
C to , qo := hmmf and C to , qo := lrmmf . 

Then Ct 0i9o = Ct 0) q . 

Proof: Let to G I and go G C(to), the first inclusion Ct , go C Ci 0)(JO is obvious and we only deal 
with the other one. So let us fix v G Ct , 9o - By definition, there exist vectors v^q G c ( t+ ^~i 



such that 



v = lim 1 3 . 



We write 



€ Pc(t +fc)(go + ^h,t,g) and tu/^g := ^ G 



h h 

It comes 

|ft^,t,g ~ Zh,t,q + qo\ _ d C(t +h){qO + ft«ft,t,g) 

h ~ h 

Moreover g + - z hM £ N(C(i + h),Zh,t,q)- Let £ A)tjg G Pc(to+/i)(<? + hv h,t, q ), we have 

1 90 + ~ ^M,?! - \€h,t,q ~ Zh,t,q\ + \Ch,t,q ~ (?0 + ft«h,t,«)l- ( 16 ) 



K,t,<? - wm, 9 I = — j: — = j: ■ ( 15 ) 



8 



FREDERIC BERNICOT & JULIETTE VENEL 



Then for h > 0, let us choose t/j G Z such that \th — to| < ^ 2 anci % £ C(th) satisfying 
\qh — 9o| < h 2 . So by this way, 

\€h,t h>qh - (90 + hv h,t h , qh )\ < \ih,t hm - {Qh + hv h)th) g h )\ + \q h - 90 1 

< d C (to+h)(qh + hv hjth ,q h ) + \qh ~ go I 

< cfo(C(t + ^Cfe + /i)) + |g/, - g 1 

< co\t h - t 1 + |g h - g 1 < (co + I)/* 2 - (17) 

In addition, 

^c(t +ft)(go + hv h>th , qh ) < d C ( to ){qo + hv h:thtqh ) + d H {C(t + h),C(to)) 

< h\v h)th , qh \ + c h 

and 

dc(t +h)(qh + hv h>thjqh ) < d c{th) (q h + hv hitht9h ) + d H (C(t + h),C(t h )) 

< h\v h>thm \ + co|th -i - 

For /i small enough, it comes 

d C (t +h){qo + hv hjthjqh ) < I and d c[tQ+h) {q h + hv h ^ qh ) < -. 

By Theorem 4.8 in [8], since C(to + h) is ?7-prox-regular, the projection operator is Lipschitz on 
a neighborhood of C(to + h). As a consequence, for h small enough 

\€h,t h ,q h ~ Zh,q h \ = \Pc(t +h)(qh + hv htthj q h ) - ~Pc(t +h)(qO + hv h ^q h )\ 

< 2r] II 

- 2r] - d c{to+h) (q + hv h , th ,q h ) - d c{ t 0+h) (q h + hv htth , qh ) 

< 2\q h - 90 1 < 2h 2 . (18) 
It follows from (16), (17) and (18) 

1 90 + h Vh,t h ,q h ~ Zh,q h \ < \£h,t h ,q h ~ z h,t h ,q h \ + \€h,t h ,q h ~ (go + hv h ^ qh )\ 

< 2h 2 + (co + l)h 2 < Ch 2 , 
for some numerical constant C. It follows from (15) that 

\Vh,t h ,q h ~ w h , th ,q h \ < Ch 

and so 

v = lim v h t q = lim w h t hAh 

with Wh,t h ,q h G Cfa+fe) »?o ^ jj ence v £ Ct 0>qo , which concludes the proof of this inclusion Ct 0tqo C 
Cto,qo- Q 

2.2. The sets of admissible velocities. — In this subsection, we consider a set-valued map 
C(-) Lipschitz continuous and taking (nonempty) uniformly prox-regular values. We aim to 
describe the sets of admissible velocity Ct <q with the help of the whole set 

tt : = {(t, q), tel and q G C(t)} 

and to prove its convexity. 

Let first recall the notion of the derivable tangent cone (see Section 6. A in [36]). 

Definition 2.10. — Let Q be a closed set ofM. d . A vector v G Tq(x) is said derivable if for all 
small enough e > there exists v e G Q=2L suc fi tfi a t Ve converges tov. We note Tg(x) C T Q (x) 
the set of derivable vectors, which can be seen as the following inner limit 

O — x { O — x 1 

TR(x) := liminf — = <u£ M. d , v = lim v e with v € G — > . 

QK ' e\o e \ e\o e J 
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A closed subset Q is said geometrically derivable if every tangent vector is derivable, i.e. for all 

x g Q 

T Q (x) = Tg(x). 

In this particular case, the inner and outer limits are equal and so the limit is well-defined: 

Tg(s)=Tg(s) = fca^. 

For example, it is well-known that every uniformly prox-regular or even Clarke's regular set is 
geometrically derivable (see Corollary 6.30 in [36]). 

These definitions allow us to describe the set of admissible velocities with the the tangent cone 
of f2. For (t, q) G f2, we have defined the set of admissible velocities 

r J r e C(t + e)- g- 

(-t a ■= i v = nm v t , with v e G 

I e\0 e 

O 

Proposition 2.11. — For every (t,q) £ f2 with t G X, we have 

C t , q = {ueR d , (l, u) £Tg((t, ?))}. (19) 

Proof: Let us denote by Ct >q the set of the right side in (19). It is obvious that for v G Ct j(? , 
the vector (l,t>) is tangent to f2 at (t,q) and derivable. Consequently the inclusion Ct <q C Ct i9 is 
proved. 

Let us now study the inverse inclusion. So consider v G Cf i9 . By definition (l,v) is a derivable 
vector to O at so for small enough e > 0, there exists a vector (s e ,v e ) converging to (1, u) 
such that 

(s e ,v e ) G 

which is equivalent to 

q + ev e G C(t + es e ). 

Since the set-valued map C is Lipschitz continuous (with a Lipschitz constant Co), it comes 

dc(t+ £ )(9 + «0 < djy(C(t + es e ),C{t + e)) < c e|s e - 1|. (20) 
So let us denote z e G Pc(t+e)( ( ? + ev e) and choose u; e such that 

z t := q + ew £ . 
By definition, z £ G C(t + e), which means 

n-(t,g) 
(l,«; e ) G ^ . 

Moreover (20) yields 

\v e — w e \ < co\s e — 1|. 

Since s e converges to 1 and v e to v, we deduce that w t converges to v too. Thus v G C^ g , which 
ends the proof of the inclusion Ct j(? C C^ q . □ 

Then we can state the main result concerning the sets of admissible velocities. 

o 

Proposition 2.12. — For every (t,q) G with t El, the set Ct, q is convex. 

Proof: If Ct^q is supposed to be nonempty, we can choose two vectors v 1 and v 2 belonging to 
Ct t q. Let a G]0, 1[, we aim to prove that 

v := av 1 + (1 — a)v 2 G Ct A - 

According to the previous lemma, for i = 1, 2 and small enough parameter e, we have a sequence 
v l t , converging to v l such that 

q + ev\ G C(t + e). 
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Writing v e := av\ + (1 — a)v 2 , we want to estimate dc(t+€){l + ev e)- Suppose that q + ev e does 
not belong to C(t + e) and let z e G P C ( t+e )((| + ev e ). So q + ev e — z e is a proximal normal vector 
at z e . Since C(t + e) is ry-prox-regular, thanks to the hypomonotonicity property of the proximal 
normal cone (see Proposition 2.6), we have for i = 1, 2 

, . „ \q + ev e — zJ . , l2 
(q + ev e - z e ,q + ev e - z e ) < \q + ev e - z e \ . 

Multiplying by a the previous inequality for i = 1 and by 1 — a the one for i = 2 and summing 
them, we obtain 

\q + ev e - z e \ < -3- \a\q + evl - z € \ 2 + (1 — a)\q + ef^ - zJ 2 l . 

Since q + ev\ € C(t + e) for z = 1, 2, we have 

|<7 + ei)* - z e | < \q + ev e - z £ \ + e|i) e - u*| 

< 2e\v e -vl\ 

< 2e(\v l e \ + \v 2 e \) < 4Me, 

where M denotes a uniform bound of v\ and v 2 (since they are convergent). So we deduce that 

8M 2 9 , . 

\q + ev £ -z e \< e 2 . (21) 



Then let us choose w f such that 



From (21), we know that 



q + ew e G C(i + e). 



em f — u; f < 



8M 2 2 



and so w e converges to v. We have proved that v = av 1 + (1 — a)v 2 is the limit of velocities w e 
satisfying 

g + ew e G C(t + e). 

That means u G Cj >(J , which shows the convexity of this cone. □ 

2.3. Admissibility. — After having introduced the notion of prox-regularity, let us now 
present the concept of admissibility. 

Definition 2.13. — A set-valued map C : X ^ M. d is said admissible on I := [0, T] if it takes 
uniformly prox-regular values (with a same constant) and if there exist 5,r,T > 0, and for all 
to G [0, T] sequences (x p ) p and (u p ) p with \u p \ = 1 and x p G C(io) such that for all s G [0,T] 
with |io — s\ < t, (B(x p ,r)) p is a bounded covering of the boundary dC(s) and 

Vp, VxG aC(s)n J B(x p ,2r), Vu € N(C(s),z), (v,u p ) > 5\v\. (22) 

For an admissible set, the "good directions" u p allow us to build inward cones. 

Lemma 2.14- — Let us consider an admissible Lipschitz set-valued map C : X 14 M rf and let 
us keep the notations of the previous definition. 

Fix to G L, x p G C(to) and v > 0. Then ift£L with \to — 1\ +u < r and x G B(x p , 3r/2) f] C(t), 

B ( x - vkou p , ^ ) C C(t + v), 



r/5 r 

as soon as v < u min := min < — — ——r , — — r } and k := ^ + 1, w/iere c is 

[ (2k + 2c + dy 2(c + + 2k ) J 

the Lipschitz constant of the set-valued map C(-) (see (8)) and r] the prox-regularity constant. 
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Proof: Let v < v m in and x £ B(x p , 3r/2) f]C(t). We consider z = x — vkqu p + 6, with \9\ < ^ . 
Suppose that z ^ C(t + v) and set y £ Pc(t+v)( z )i v := z — y £ N(C(i + v), y) \ {0}. Necessarily, 
y belongs to B(x,r/2) because 

\x — y\ < \x — z\ + \z — y\ < 2\x — z\ + cqv < 2 ( kqv + — J + cqv < v(5 + 2kq + cq) < r/2, 



where we have used that d c ^ t+u ^ (z) < dc(t)( z ) + d>H(C(t),C(t + f)) and x £ C(t). 
Letf; G P(;(( +[ ,)(i), since C(i + z^) is ry-prox-regular, it comes 

(£ -y,v) < - y\ 2 \v\. 



That yields 



277' 

277 
1 



(£ - z>«) + M 2 < ^- (|£ - »| + \ x - y\f \ v \ 



< ^- v 2 {5 + 2k + 2c ) 2 \v\. 
lr\ 

Consequently, 

k (vu p ,v) - (6,v) 

1 _ r 1 .1 

■ (23) 



< + 2k ° + 2c o)) 2 M + d C (t+u)(x)\v\ < 



^-(v(5 + 2k + 2c )) 2 + c v 
2rj 



As y £ B(x,r/2), y belongs to B(x p ,2r) and the admissibility assumption gives (u p ,v) > 5\v\. 
Thus (23) implies that 

k Sv - v^- < —{y{8 + 2k + 2c )) 2 + c v 
2 2r/ 

and so that v > (2k +<5+2c )^ ' wn i c h leads to a contradiction. Thus z £ C(t + v). □ 
The following lemma is a consequence of the previous one. 

Lemma 2.15. — Under the previous assumptions and notations, fix to £ I and x p £ C(to). 

C(t + v) x 

For all t £ X and x £ B(x p , 3r/2) f~| C(t), for all v < v m in, ~ K oU p £ . 

Moreover for all y £ c(t+ ^~ x f] B(0, £), 

-k u p £ T C (t W -, (y). 

In addition, —kqu p — y belongs to Ta(t+v)-x (y) as soon as 

v\y\ 2 < &q/2. (24) 
Proof: Concerning the first point, by Lemma 2.14, x—kqvu p £ C(t+v) if x £ B(x p , 3r/2) f] C(t) 
and v < iSmin- Thus —kqu p £ ^— ^ — — . Let y £ c '( t+ ^~ x p B(0, ^), x + vy £ B(x, r/2) C 

B(x p , 2r) and x + vy £ C(t + v). The admissibility assumption implies that for all v £ N(C(t + 
v),x + vy), 

(—Up, v) < —5\v\ < 0. 
Since N(C(i + v), x + vy) = N ( c(t + t/) " x ,y), we have for all v £ N( C(f+ ^~ x , y), 

(-u p ,v) < 0. 

As a consequence, — u p £ ~N( C( - t+ ^~ x ,y)° = T c(t+ v )-x (y) (due to Lemma 2.8). 

Let us now prove the last point. Since c ( t+ ^ x is 2_p rox _ re g U i ar anc j contains a £ Pc(t+1 ^^ 
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(satisfying \a\ < c ), for all y G C(t + v) and u G N( c{t+ ^~ x ,y) = N(C(t + u),x + uy) 

(a - y,v) < ^-\v\\a - y\ 2 < -\v\ (c§ + \y\ 2 ) , 
zrj r] 

by the hypomonotonicity property of the proximal normal cone. It follows from (24) that 

(-K u p -y,v) < -k S\v\ + (-y,v) < -k S\v\ + -\v\ (cq + \y\ 2 ) + c \v\ 

3(5. v. l9 

< \v\ + - \v \\y\ z < 

4 rj 

(since vc\ <r)6/4 and kq5 = 5 + cq), which proves the expected result (thanks to Lemma 2. 



□ 



Proposition 2.16. — Let C(-) be an admissible and Lipschitz continuous set-valued map on 
1. For all t$ £ I and qo G C(to)> we se t 

r .. . - C(t + h)-q ~ C (t + h)-q 

w , qo := limmf and tt (hqo := limml 



h\o h h 



?ec(t) 



Then C tom = C tom / 0. 

Proof: We refer the reader to Proposition 2.9 for the equality Ct 0t q = Ct 0) q - It remains us to 
check that these sets are nonempty. 

If qo G Int(C(to)) then we easily have Ct 0j q = Ct 0j q = M d . Else go S dC(to) and by Lemma 2.14, 
there exist u p ^ 0, > and v m i n > such that for all v g]0, v m in[, Qo — vkqu p G C{to + v). 
So we deduce that — kqu p G Ct , qo / 0. □ 



3. Discretization and convergence of approximate solutions 

This section is devoted to the proof of Theorem 1.3. As usual, we obtain existence results for 
(9) by proving the convergence of a sequence of discretized solutions. 

To do so, we extend the Catching-up algorithm (proposed by J.J. Moreau for the firt-order 
differential inclusions) to the considered second-order problem. Let us describe the numerical 
scheme. 

Let h := T/N < Vminl^ be the time step, where v m i n is defined in Lemma 2.14. We denote by 
qj^ G M d and G M d the approximated solution and velocity at time £JJ = nh for n G {0, .., N}. 

The approximated solutions are built using the following scheme: 

1. Initialization : 

(ql ql) ■= (qo, qo + huo + h 2 /°) with f h := i J** f(s, q° h )ds. (25) 

Since qo G Int(C(0)), q\ belongs to C(h) for h small enough (such that (\uo\ + j x F(t)dt)h < 
ddC(o)(qo)) ■ 

,n + l 
1 f*h 

2. Time iterations: q l h are given for i G {0, We define fh := T / fi s ,qh)ds and 

h 

1h +1 e P c( ^i) [H - qV + h 3 ffi . (26) 

This algorithm is a "prediction-correction algorithm": the predicted point 2q^ — q^~ l + h 2 f%, 
that may not be admissible at the time t/J +1 , is projected onto C(^ +1 ). 
We define the following functions: for all t G 

n n+1 - n n 

q h (t) := ql + (t - t n h ) Qh - Qh (27) 
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and the velocity Uh by 



n+l _ n 

u h (t):=u n h +l := qh - qh . (28) 



Note that, for every h < ^^^p^dt ' ^ ne sc h eme is well-defined. Moreover the computed 
configurations are feasible : 

V/ i >0,VnG{0,...,iV}, q h {tl) = <ft G C{t n h ). (29) 

For the intermediate times t G]i^,i^ +1 [, the point qh(t) may not belong to C(t). However from 
(29) and the Lipschitz regularity of the set-valued map C(-) (see (8)), we have the following 
estimate : 

Vh > 0,Vt G X, d c(t) (q h (t)) < m a x{d c{t) (q%),d c{t) (qZ +1 )} < c h. (30) 

The proof of Theorem 1.3 is quite technical and will be decomposed into 4 steps, which we shall 
briefly describe below. 

— In Subsection 3.1, we obtain uniform bounds on the computed velocities Uh when h goes 
to in L°°(X,R d ) and in BV(I,R d ). 

— In Subsection 3.2, we use compactness arguments in order to extract a subsequence of 
(qh,Uh)h>o converging to (q,u). 

— In Subsection 3.3, we check that the limit functions q, u satisfy the second order differential 
equation (the momentum balance) and the initial conditions. 

— In Subsection 3.4, we verify the impact law for u. 

After that, we will have proved that the limit function q is a solution of Problem (9), which 
shows Theorem 1.3. □ 

3.1. Uniform estimates on the computed velocities. — This subsection is devoted to 
the proof of uniform estimates in BV(X, M. d ) for the computed velocities in order to extract a 
convergent subsequence by compactness arguments. 

First, we check that the velocities are uniformly bounded in L°°(I, M. d ). Aiming that, we prove 
the following lemma giving a first estimate on the velocities. 

Lemma 3.1. — For all integer n G {0, ...,T/h} 

K +1 |<2K + ^| +C0 . 

Proof: By rewriting (26) in terms of velocity, we deduce that 

< +1 e P C (^)-, g K + hfn . 



With z := P c(t n+i } (g£) (as c h < co^mm/4 < rj, the projection de q^ is single- valued), it comes 

K + ^4 n -< +1 l< K + hfj:- 



h 

The proof is also achieved thanks to \z - g£| < d H (C{tl +1 ), C(tD) < c h. □ 

To iterate this reasoning, it would be very interesting to obtain a similar estimate without the 
factor 2 in Lemma 3.1. This is the main difficulty in order to obtain a uniform bound of the 
velocities and we solve it in the following proof. 

Proposition 3.2. — There exists hi > such that the sequence of computed velocities (uh)h<h! 
is bounded in L°°(X, R d ). We set 

K := sup ||u/i||i°°(x) < °°- 

h<hx 
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Proof: 1 — ) Estimate on the velocities for small time intervals. 

Let us fix h < hi with h\ later defined and satisfying h\ < min{z/ m j n , r}. Consider a small time 
interval C I of length satisfying 



h<\t+-t-\<mm{ — ^,-} (31) 




where kq is introduced in Lemma 2.14 and uq is the smallest integer n such that t 1 ^ > i_. Thus 
G We are looking for a bound on the velocity on this time interval. Rewriting the 

scheme in terms of velocities, we have that 

h 

involving u n h + hf% - u n h +1 G N ^^"^ , u™ +1 

By the admissibility of C, there are balls B(x p ,r), which cover C(t^°). So there exists at least 
one index p such that q^° G B(x p , r) and we denote uj := u p . 

Then for n satisfying tjj G t+] (since h+\t r J l — < r), Lemma 2.15 yields that for h < v m in 
if 

q n h eB{ql\r/2) (32) 

(which implies q^ G -B(x p , 3r/2)) and 

h\u n h +1 \ 2 < V 5/2, (33) 

then -k w - < +1 G N , v%+A . So we deduce that 

{u n + hf n_ u n+l^_ KoW _ u n+l ) < 0j 

which implies 

K +1 + k w\ < \ul + hf% + K ^| < K + K o w \ + h \fh\- 
We set m the smallest integer (bigger than no) such that m + 1 does not satisfy (32), (33) or 
^m+i ^ By summing these previous inequalities from n = tlq to n = p with uq < p < m, 

we get 

Vp G {no, m}, + K oH < + K o^| + / F(t)dt. 

Jo 

Finally, it comes 

/■C' +1 

sup \u p h \ < \ul°\ +2k + / (34) 
:p<m+i Jo 



no<p<m+l 

By integrating in time, we deduce 

^m+l _no I ^ | |„,no 



< ( |tt£°| + 2K + jf " j - t_| + /l) < y < 



2r r 
2 : 



by (31). Consequently by definition of m, t™ < t+ < t™ +l as soon as (33) holds for n = m + 1. 
Moreover thanks to Lemma 3.1, we have 

K +2 \ < 2|n™ +1 | + c + 2 / F(t)dt. 

Jo 

From (34), (33) is satisfied for n = m + 1 as soon as 

h^2{\u n h °\ + 2n Q + £ F(t)dt)+c + 2 £ ' F(t)dtj < ^. 
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Finally, we obtain 



sup K| < |u£°| + 2k + / F(t)tft, 
<t£<t+ Jo 



as soon as 

h \2(\ul 1 + 2k + J^ F{t)dt)+c + 2 J F(t)dt\ < ^. 
2—) End of the proof. 

Let h < hi (later denned), we are now looking for a bound on Uh on the whole time interval 
I = [0, T], in assuming that T < r/2 without loss of generality. Let us start with t- = t(0) := 0. 
We know that with 

r 



where 



we have 

0<t£<t(l) 



t + =t(l) :=min{— ,T 



A(l) := \u \ + 2k + [ F{t)dt 
Jo 

sup K| < A(l) < \u Q \ + 2k + f F(t)dt, 
t?<t(D Jo 



as soon as 



h(2A(l)+co + 2 jT F(t)dt\ < ^. 



Then, let us suppose that there exists n x such that t(0) < t^ 1 < t(l) < t™ 1+1 . We have 
< Si := t(l) — t^ 1 < h. In that case, we set t- = i^ 1 and 



t + = t(2) := min <^ t? 1 + 



5,4(2)' 



min^(l)-5 1 + ^,T 



with 



A{2) := |uq| + 4k + 2 f F(t)dt > M^ 1 )] + 2k + f F(t)dt. 

Jo Jo 
From the previous point, we deduce that 



SUp \Ufr\ < SUp \Ufr\ 

t{i)<q<t(2) C<q<t(2) 

rT 



as soon as 



Hence 



<\ul 1 \+2n + [ F(t)dt 
Jo 

< \u \ + Akq + 2 / F(t)dt = A{2) 
Jo 

h Ua{2) + c + 2 J F(t)dt\ <^. 

sup \ul\ < A(2) = \uq\ +4k + 2 f F(t)dt. 
:t"<t(2) Jo 



0<*™<t(2) 

By iterating this reasoning, for any integer k > 1 we set 

rT 



A(k) := \uq\ + 2kn + k [ F(t)dt = A(k - 1) + 2k + f F(t)dt 

Jo Jo 
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and 



t(k) := min <^ t(k - 1) - 5 fc -i + 



5A(k) 



k-1 



min 



i=l 



where 5k < h for all fc. This construction of can be made while t(k — 1) — i(/c — 2)>h. This 
condition will be verified as long as 

r 



Sk-2 + 



> h. 



5A(k - 1) 

Using the fact that < Sk-2 < h, we see that we can construct t(k) for k < N verifying 

r 

> 2h, 



5 ( | i*o | + 2(& - 1)k + (k - 1) J T F(t)dt 



which is implied by 



k < k (h) :-- 



1 + 



-i 



(35) 



Consequently, we know that the velocities can be bounded on [0, t(ko(h))] as follows 

sup \u k h \ < A(k (h)) 
o<ti<t(k {h)) 



where 

under the property 
/; 



feo(»-l k (h) 



t(k (h)) = min I - 6 i+Yl 



1 5A(i) 



(36) 



F(t)dt < 7]5/2. 



2 sup \ul\+c + 2f F(t)dt\ <h(2A(k (h)) + c + 2 [ . 

0<t™<t(k (h)) Jo J V JO 

Now, using the fact that ko(h) goes to infinity when h goes to zero, that the harmonic serie 
diverges and that (35) yields 

k (h)-i 

Y, ^ <hk (h)<c, 



■i=i 



we see by (36) that t(ko(h)) is equal to T for h small enough. Therefore, there exists ho < v m i n 
such that, for h < ho, T = t (Icq (ho)) = t(ko(h)). Finally, we see that, for h < ho, t(k) can 
be constructed until k = ko(ho) (since h — > ko(h) is non decreasing) and that t(ko(ho)) = T 
(independently from h < ho). Hence Uh can be bounded as follows 



sup \ul\ < A(k (h )), 



if 



h ^2A(ko(ho)) + c + 2 F(t)dt\ <i]5/2. 



(37) 



(38) 



However, (38) holds for 

h < hi := min < 



ho, 



T]6 



2(2A(ko(ho)) + co + 2f Q T F(t)dt 
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So we finally obtain the uniform bound 

sup sup {u^l < A(k (h )) = \u \ + 2k (h )K + k (h ) / F(t)dt, 
h<hi 0<t£<T 70 

which concludes the proof of a uniform bound in L°° for the velocities u^. □ 

Having obtain a uniform bound of the velocities, we can now prove that they have a uniformly 
bounded variation on the whole time-interval /. 

Theorem 3.3. — There exists h,2 G]0, h\ \ such that the sequence of computed velocities (uh)h<h 2 
is bounded in BV(I,R d ). 

Proof: We adapt the proof of Theorem 3.2 in [5]. To study the variation of Uh on /, we split I 
into smallest intervals. We define (sj)j for j from to P such that: 

s = , s P = T, 

\s j+ i - Sj\ = -minjr, -^j , for j = . . . P - 2, 



\s P - sp-i\ < - min |t, , 



where r and r are given by Definition 2.13 of the admissibility and K is the Lipschitz constant 
of qh (see Proposition 3.2). All these constants are independent on h and such a construction 
gives 

2T 



P 



mm 



+ 1, (39) 



which is independent of h. Then, for all h < h±, we define n 3 h for j from to P — 1 as the first 
time step strictly greater than s,-: 



,n J h — 1 - n J h 
*ft — S 3 < *ft ' 



and is set equal to N (t^ = t h h = T). 



In the sequel, we suppose h < min{|sj + i — Sj\}/2. We also obtain a strictly increasing sequence 

j 

of (t h h )j verifying 



lt?-^ rl |<min{r,-^}. (40) 



The variation of on I can be estimated as follows 

N-l P-l 
Var/(u fc ) = K +1 ~Uh\ = Yl Var i u h 

n=0 j=0 

where 

n{ +1 -l 
Va rj K):= ^ 

n = n h 

corresponds to the variation on [t h h , t h h [. To study these terms, we recall that 

.n+l 



< +1 eP C(i » + V, r [< + hm (41) 



by writing the scheme in term of velocities and we state the following lemma: 
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Lemma 3.4- — There exist uniformly bounded vectors y n h (\y n h\ < kq) such that, for all small 
enough h, for all j G {0, . . . , P} and n G N H [n^, [ , we have 



xi €P C{t ^+i ) _ q n[x ] =>- \xi - x \ < I (\x - y n h\ 2 - \xi - y n h\ 2 



as soon as x\ is bounded by K (constant introduced in Proposition 3.2). 
Proof: First for small enough h < v m in and h < r/2, it comes 

n G [n{,n{ +1 [=^ \tf - t n h \ + h < r/2 + r/2 = r and \q< - tf\ < K\t< - t n h \ < r/2 

7 7 

So thanks to Lemma 2.14 (applied for t h h ,q h h ), we know that there exist unit vectors v n h such 
that 

ne[nln{ +1 [ B(-k v<,5/2) c U - ' Qh . (42) 



Indeed v n h is "a good direction" , given by the admissibility assumption, associated to the point 

j 3 

(V>9h fc )- 

Then, we develop similar arguments than those used in [10] and [5, 31, 13]. In these works, 

C(t" +1 ) — q n 

the set onto which the velocity is projected was convex. In the present case, the set h h - 

is ?///i-prox-regular, which is slightly weaker. 

j 74 



Let n belong to [^,n^ +1 [. We define 



2: n h : = y n h + where y * := -kqv 

2 \xq — xi\ 

(Here we suppose xq 7^ x±, else the desired result is obvious.) We have 

C(t n h +1 )-g n h 



z n * G B(-K Q v n h,5/2) C 



h 



The point xi being the projection of xq onto the ^-prox-regular set — h -j i -, we get 

(x - xi,z n h - xi) < ^-\x - xi\\z n h - xi| 2 
lr\ 

thanks to the hypomonotonicity property (see Proposition 2.6). Thus 

\x - V n ' h ? = \xi - y n h\ 2 + \xo - xi| 2 + 2{z< - i/"Mo - x x ) + 2(a?i - z<,x Q - xi) 

> \x\ - y n h\ 2 + 2{z n h - y^^XQ - x\) - -\x - xi\\z n ^ - xA 2 

V 

1 1 2 rl 1^1 11 1 2 

> \xi — y h\ + d\XQ — xi\ \xq — x\\\z h — x\\ . 

V 

Using that the vectors z n h are uniformly bounded by kq + 5/2 and that x\ is bounded by K, it 
follows that for h < 2{rq J/ 2+k) -2 



ko-y n H > |xi - y n '^ + -|x - 

This, together with the fact that the vectors y n h are uniformly bounded by kq, ends the proof 
of Lemma 3.4. □ 
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We now come back to the proof of Theorem 3.3. For n in [n? h , using (41) and the previous 

lemma (with xq = u\ + hf£ and x\ = u% +1 ), it comes 

2 



u n h +x -u n h -hf%\ < - 5 [\x G -y<\ 2 -\ Xl -y<\ 2 ) 



6 

< l(K + hfZ-y n i\ 2 -\ul +1 -y<\ 2 

< - 5 {K - y<\ 2 - K +1 - y<\ 2 ) + - g \hfZ\ 2 + - 5 \hm\u n h - y<\ 

< | (K - y ni \ 2 - K +1 - y ni \ 2 ) + - § \hm 2 + - 5 \hm (k + «o) . 



By summing up these terms for n from to — 1 we get 

Var,M= £ |„^-<| < ^(|u?-^| 2 -|^ +1 -^| 2 )+ £ | W | 



and finally 

P-l „ P-i 



Var(^) = ^Var J K) < ^ (V? - y<| 2 - - y<| 2> ) + ^ (7 F(t)dt) 

j=0 j=o ^ ' ' 

This completes the proof of Theorem 3.3, since P does not depend on h from (39). □ 

3.2. Extraction of a convergent subsequence. — The previous uniform bounds on the 
computed velocities allow us to extract a convergent subsequence. 

Proposition 3.5. — There exist subsequences of (q^) and (u^) (still denoted by qh and Uh) 
which respectively converge to q £ W l '°°{T,M, d ) and u £ BV(l,R d ). Moreover (u h ) h strongly 
converges to u in L l (T, M. d ) and for all t € X 

q(t) € C(t). (43) 

The initial condition is satisfied: q(0) = qo and u(0) = uq. 



Proof: By Proposition 3.2, the sequence (uh)h is bounded in L°°(I, Mr). Arzela-Ascoli theorem 
proves that (qh)h is also relatively compact in W 1,OQ (l, M d ). So up to a subsequence, we can 
assume that qh strongly converges to q G W 1 ' OD (I,'R d ). 

Moreover, as (it^) is bounded in BV(Z, M d ), there exists a subsequence (still denoted by Uh) 
converging to u in L\l,R d ). It is easy to show that necessarily u = q in the distributional 
sense. In addition, by the uniform bound of the variation Var(u^), u belongs to BV(I, M d ) with 
Var(u) < sup h Vax(uh). 

Inclusion (43) is a direct consequence of the uniform convergence of qh to q together with (30). 
Let us now check the last point concerning the initial condition. Since for all h > 0, Q7i(0) = qo 
so q(0) = qo. Moreover, the initial point qo £ Int(C(0)). The maps qh and q are Lipschitz with 
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the same constant, which implies there exist s > and I > such that for all t G [0, s] and all 
small enough h > 

ddC(t)(qh(t)) + d 9c[t) (q(t)) > I. 
So the computed points are far away from the boundary of C(-) during the whole interval [0, s]. 
That implies 

Q n h +l = h\ - ?r 1 + * 2 #. 

in other words 

< +1 = tfl + 
for n verifying G [0, s]. So we deduce that 

K +1 -ul\< / F(t)dt, 
which gives by summing all these inequalities 

rs+h 

sup \uh(t) — uq\ < / F(t)dt. 
te[o,s] Jo 

Since Uh converges almost everywhere to u (up to a subsequence), it comes that 



\\u - uo\\ loo( t s] K d) < / F(t)dt. 

Jo 

Taking the limit when s goes to implies the desired result: it(0) = Uo. □ 

3.3. Solution of the continuous differential inclusion. — In this subsection, we prove 
that the limit function (obtained in the previous subsection) satisfies the differential inclusion 
of Problem (9), according to Definition 1.2. 

Proposition 3. 6. — The limit function q satisfies the continuous differential inclusion: there 
exists k G BV(I,M. d ) such that in the sense of time-measures 

dq + dk = f{t, q)dt (44) 

and the differential measure dk is supported on {t, q(t) G dC(t)} : 

\k\{t) = I \{s)edC{s)d\k\{s), k(t)= [ Z(s)d\k\(s), (45) 

J J 

withC(s) G N(C(s),q(s)), |£(s)| = 1 and \k\(t) := Var(k, [0, t\). 
The idea is to let h goes to in the relation 

< + hfZ - ul + 1 g N (46) 

We refer the reader to [20], [37] and [6] for similar reasonings in the framework of first order 
differential inclusions. 
Proof: The scheme 



< +1 eP c(t „ +v „K + ^] 



implies (46) and so 



h 

< + hft - u n h +1 G N(C(^ +1 ), q n h +1 ). (47) 



Let us define a piecewise-constant function kh, defined for t G [t^ 1 ,t^ l [ by 

n— l 



k h (t) :=k%:=u -ul + y2 1 ^h)ds. 

i=0 17 t h 

So we have 

- K = < + hfft - u n h +1 G N(C(i£ +1 ), <£ +1 ). (48) 
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By Proposition 3.5, converges to u in L 1 (X, M. d ) and qh uniformly converges to q. Moreover 
/ is Lipschitz with respect to the second variable thus we deduce that kh strongly converges to 
some function k £ L l {X, M. d ) verifying for almost every t £ X 

k(t) := u - u(t) + / f(s,q(s))ds. 



Thanks to the uniform bounded variation of Uh, k belongs to BV(T, M, d ) and then (44) holds. 
So it remains us to check (45). We recall that f2 is the set of (t, q) with t £ X and q £ C(t). Let 
X ■ X x M. d — > M. d be any nonnegative smooth function compactly supported in Int(f2). Then 

0< / X (s,q(s))d\k\(s) <liminf / X ( a , 9ft ( 8 ))d|fc ft |(*). (49) 
Jo Jo 

However, 

T N-i q+i 

/ x(-,flfc(-))d|*fc|(*)=E / X{s,q n h +{s-tl)ul +1 ){kl +l -k n h )ds. 

Assume that for some integer n, k^ +1 — k% ^ 0, then + + h 2 f% g C(i^ +1 ). Thus, for 
every i£ 



d aC (t)(Qh)<h(co+K + J F(*)di 



because q% £ C(iJJ), \u^\ < K (see Proposition 3.2) and Co is the Lipschitz constant of C. 
Consequently, for h small enough 

ddc(t)(Qh) + hK < d H (supp( x ), 

and so q% + (s — tfyu 1 ^ 1 ) = 0, for s £ [iJJ, Finally, we obtain that for h small enough 

/ X(s,qh(s))d\k h \(s) = 0, 
J o 

hence from (49), we obtain 

X (s,q(s))d\k\(s)=0, 



for any nonnegative smooth function Xi compactly supported in Int(Q). Taking a sequence of 
such functions converging (increasingly) to li n t(n), then we have 

hnt(n)(s,q(s))d\k\{s) = 0, 

o 

or equivalently 

\k\ T = [ lan(s,q(s))d\k\(s)= I l dc(s) (q{s))d\k\(s). 
Jo Jo 
To finish, it remains us to check that ll dk(s) £ N(C(s), q(s))" . We move the reader to [20] 
(the end of the proof for Theorem 1.1 in [20]) for precise details. Indeed the arguments relie 
on the hypomonotonicity property of the proximal normal cones (see Proposition 2.6). By (47), 
u h "t" hfh ~ De l° n g s to N(C(i^ +1 ), q^ +1 ) for all integer n. So for every continuous map 

<j) : X ->■ R d such that <f>(t) £ C(t), we have 

mi +l ) ~ <f h + \< + hfjt - u n h +l ) < + fctf - < +1 ||^(^ +1 ) - ^ +1 | 2 . 

Summing all these inequalities from to N — 1, it comes with (48) for every nonnegative function 



1 



T 



[ m^(0h(t)) - q h (e h (t)),dk h (t)) < — I m\Hdh(t)) - q h {6 h (t))\ z d\k h \{t) 
Jo 2 V Jo 
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where we denote 6 h {t) = for t G [t£,tE +1 [. Since dkh and d\kh\ are uniformly bounded 

measures. Up to extract a subsequence, we can assume that they are weakly convergent to dk 
(the differential measure of k) and da (where da is a nonnegative measure). Necessarily, the 
measure dk is absolutely continuous with respect to the measure da. So there exists a bounded 
and measurable function g such that dk = gda. Taking the limit in (50) when h goes to 0, it 
comes 

T I rT 



mm) - q(t),g(t))da(t) < — / mm) - <K*)r<M*). 

2?? Jo 

Since this inequality holds for every nonnegative function ip, we deduce that for all t G X and 
all map (ft 

{<i) (t)- q (t),g(t))<^(t)- q (t)\\ 

which yields by Proposition 2.7 that g(t) G N(C(t),q(t)). Indeed for every to £ I and <fro G C(to), 
there exists a continuous map 4> : X — > M. d satisfying 

f <f>{t) eC(t), Vtel 

[ <f>(to) = <?V 

It suffices to consider the solution of the following sweeping process (see [41]) 

J -0(t) G N(C(t), 0(t)) for a.e. t G X 
i <K*o) = 00- 

That also concludes the proof of (45). □ 



3.4. Collision law. — Finally, Theorem 1.3 will be proved, provided that the collision law is 
satisfied for the limits u and q, which is the aim of the following proposition. 

Proposition 3.7. — The impact law is satisfied: 

Vt GX, u + (t )=P CtoMto) (u-(to)). 

Proof: Note that, from Proposition 3.5, u G BV(I, M. d ), so that the left-sided u~(to) and the 
right-sided u + (to) limits are well-defined. 

The proof is quite technical so for an easy reference, we remember the definitions of the sets 
Ct 0t q (see Definition 1.1) 

C *o,9o ■= <v = lim v e , with v e G 

Moreover, we recall that these sets are nonempty due to Proposition 2.16. From now on, let us 
fix the instant to G X. The desired property 

n + (io) = P Cto , 9(to) (^(i )) (51) 
can be seen as the limit (for h going to 0) of the "discretized property" 

«; +1 £p < K Wh+hr\. (52) 

h 

Step 1: We claim that 

u+(t ) G C toMtQ) . (53) 
By definition, q G W^°°{1, R d ) and u G BV(l,R d ) so we have that 

+ ,. \ r u , \ v g(^Q + g) - g(t ) 
u T (to) := hm u(to + el = hm . 

e->0+ e->0+ £ 
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The last equality comes from 

j(t + a) - q(t ) 



u+(t ) 



Since q(to + e) G C(to + e) 



«+(t ) 



u(s)ds 



to 



< sup |it + (to) — u i s )\ 

s£[to,to+a] 



-> 0. 



n + (£o) G nm 



C(t + e)-q(t ) 



C 



*0><K*o)' 



which completes the proof of (53). 



Let us now come back to the proof of (51). As we just proved u + (to) G C to , g (t ) and since C t0)9 (t ) 
is a convex set (see Proposition 2.12), (51) is equivalent to 

Vu; G C t0)?(to) , (u-(to) - u+(i ),™ -u+(t )) < 0. (54) 
So, in the following, let us fix w G C to , g (i )- In Step 2, we construct a family of points for 
v > such that w„ tends to i« when v goes to zero satisfying w„ G c ( t +^) g f or /j sufficiently 
small and (t, g) close to (io, q(to))- Then in Step 3, for each is, we go to the limit on h, t and q 
to show that (u~(to) — u + (to),w u — u + (to)) < and finally, we make v go to zero to conclude. 

Step 2: From the admissibility assumption, there exist a neighborhood U C SI C Xx R rf around 
(to,q(to)) and a "good direction" £ G R d such that for all (t, q) G U 

V«€N(C(*),g), (C,v) <-8\v\, (55) 

with a numerical constant <5 > 0. For v g]0, 1[, we consider the point w v := w + ^C- We claim 
that for each fixed v > 0, there are ej, and /i^ such that for every h < h u , (i, q) G t7, we have 

C(t + h)-q 



\t-h\ + \q-q{t )\ <e h 



w u G 



(56) 



Let us detail this point. 

Thanks to Proposition 2.16, we know that for any 9 there exists a neighborhood Vg C U of 
(to,q(to)) and /i# such that for all (t,q) G and /i G]0, hg[ 

d c(t+h)- q (w) < 9. 



Let us denote Wtq a point of P 



C(t + h)-q 
h 



(to). 



li h < h$ < 2r]/ (co + |w| + 1), it can be proved that for (i, g) G Vg, w„ G c ( t+ ^ q i. e g + hw u G 
C(t + h). By Proposition 2.7 and as d C ( t+h )(q + /uiv) < 2^, it suffices to show that for all 
x G C(t + /i) and u G N(C(i + h),x) 



as soon as 



(<? + hw u — x,v) < —\q + hw u - x\ 2 
\x — q + hw u \ < h(co + |iu| + 1). 



(57) 
(58) 



Let x G C(t + h) satisfying (58), we have 

(q + hw u — x, v ) = (q + hw + vhQ — x, v) 

= (q + /u*) ti(? - x, v) + z^(C> + h(w - w t , q , v) 

< — \v\\q + hwt q — x\ 2 — uh5\v\ + h\w — wt n\\v\ 
zrj 



1 , 

— \q + nw t , q - x\ 
2r) 



uh5 + h9 



24 



FREDERIC BERNICOT & JULIETTE VENEL 



where we have used that q + hu>t, q and x belong to C(t + h) with the hypomonotonicity property 
of the proximal normal cone and (55). Then taking 9 < min{i^5/2,ci} with c\ := l/(6(co + 
2\w\ + 2)), we get 

\q + hwt s q — x\ 2 <\q + hw u — x\ 2 + h 2 \w u — u>t,q\ 2 + 2h\w u — Wt,q\\q — x\ 



< \q + hw„ — x\ 2 + h 2 



and so 



2rj 



\q + hw 



t,q 



< —\q + hw*q - x 
2rj 

1 , 

< — \q + hw u — x\ 
2rj 

1 . 

< —\q + hw v - x\ 
2rj 



l -h 


u6 




T 


- h 


'uS 




T 



2r] 



as soon as h < v8rj. Thus, for v < c±, < min{i/5/2, ci}, h < h u := min{/ifii, v5rf\ and (t, q) G V$, 
it comes 

— \q + hwt n — x\ 2 — vh5 + h6 < — \q + hw u — x\ 2 
2t] 2t] 

for all x G C(t + /i) satisfying (58) and v G N(C(i + /i), x). That proves (57) and so (56). 
Step 3: Let us now fix the parameter v < c\. 

Thanks to the uniform Lipschitz regularity of the maps qh and their uniform convergence towards 
q, there exists h v < h v such that for e < e v /(2 + 2K) and h < h u , 

t k h ,t k h +1 G [t - Mo + e] \t k h +1 ~ t \ + \q k h - q(t Q )\ < e u . 

We recall that 

K := sup ||«h||L«,(x,]Rd) < oo. 
Consequently, as q k G C(t^), the property (56) (with i = t k ) gives wv G . ^ gh . 



Moreover 



(52) describes that 



u£ + /i/*-u* +1 €N 



C(4 +1 ) 



^ „.fc+i 



Therefore, — ^-r being ^-prox-regular, we have (due to the hypomonotonicity property) 



(4 + hf k h - u k h + \w v - u k+1 ) <±\< + hf k h - < +1 |K - u 



k+l\2 
h I • 



(59) 



We sum up these inequalities for k from n to p, integers chosen such that is the first time 
step in [to — e,to — e + h] and t p h the last one in [to + e — h, to + e]. First, we know that 



u 



k+l\ 
h l 



k=n 



rto+e+h 

< (\w\+K + l) / F(t)dt 

J tn— e 



with K := sup^ ||^||oo- We also have 

p 

S(«*-«fc +1 .^> 

k=n 



{u h {t n h -')-u h {tl),w u ). 



We deal with the remainder as follows: 

- -uj* 1 ) = j2(u k h - u k +\u k ) - \ut\ 2 + \u 



«|2 _l U,P+1|2 
h I ' 



k=n 



(60) 



(61) 
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which gives 



£<«* - 4 + \ -< +1 ) = \ E \< - < +1 \ 2 + \ [-Mr 1 )! 2 + KO 2 ] 

= IvaraK)^ ^ + ^ [-MC 1 )! 2 + MOT , (62) 

where we wrote Var2 for the L 2 -variation of a function. Using (59), (60), (61) and (62), we 
finally get : 

ivarafo)*., tP] + J [-M^r 1 )! 2 + \u h (t{)\ 2 ] + (tiftCtg -1 ) - u h {1? h ), Wv ) 



< (| W | +K + 1) F(t)dt + -J2\u k h + hf k h 

Jtn-e Zr ? 



«{ +1 lK-«s +1 i a . 



However 

E 14 + ^4 fc - < X>* - + E ^l/hl ^ Var K) + / Ti w < s i 

fe=n fc=ra fc=n ^ 

for some numerical constant, due to Theorem 3.3 and 

\w u - u k h +l \ 2 < (\w\ + 1 + K) 2 = B 2 , 
due to Proposition 3.2. Consequently, we deduce that 

ivarsK)^^ + \ [-Mtr 1 )! 2 + \u h (t{)\ 2 ] + (lifcC^- 1 ) - u fc (#,«fc) 

fto+6+h 



pto+e+n Zj 

< £ 2 / F{t)dt + —B l B 2 . (63) 



Let us now choose a sequence of e m going to zero, such that u/j pointwisely converges to u at 
the instants io ~~ e m an d to + e m (which is possible as Uh converges almost everywhere towards 
u). For each e m and h < h u , we have shown that inequality (63) holds. Then, passing to the 
limit for h — > we get 

2 Var 2(^)[t - em ,t 0+em ] + 2 H u (*° ~ €m ^ 2 + ' n ^° + £m )' 2 ] 

+ (u(t - e m ) - u(i + e m ), «;„) < B 2 F(t)dt, 

J tQ— Cm 

which gives for e m — )• 

ivar 2 (n)2_ ifo+] + ^ [-|n-(t )| 2 + k + (to)| 2 ] + («~(*o) - ^ + (*o), v) v ) < 0. 
Finally we obtain 

- \u + (t ) - u~{t )\ 2 + - [-\u-(t )\ 2 + \u + (t )\ 2 ] + (tr(to) - u + (io),Ti>„) < 0. 
By expanding the square quantities, we obtain for all v < c\ 

(n-(t )-u + (to),^-u + (t )) <0. (64) 
Recall that w v = w + v(,-> we obtain (54) by letting v go to in (64). □ 
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4. A particular case 

As explained in the introduction, all the second-order differential inclusions, already studied 
in the literature, concern a particular case where the moving set C(-) is given by a finite number 
of constraints. This section is also devoted to prove that the previous abstract result covers this 
case, as soon as the constraints satisfy some reasonable assumptions. 

So we consider the Euclidean space M. d , B the closed unit ball in M. d+1 and X = [0, T] a bounded 
time-interval. For i £ {1, ...,p} let gi : I x R rf — > R be functions (which can be thought as 
"constraints"). For t £ X, we introduce the sets 

Qi(t) := {iel d , gi (t,x)>0}, 

and the following one 

p 

Q(t) :=f]Qi(t), 
i=l 

which represents the set of "feasible configurations x" . We remember that Q := {(£, x), t £ 
X, x £ Q(t)} and we similarly define Oj in replacing Q(t) by Qi(t). Moreover we assume that 
there exist a, (3, M, k > such that gi £ C 2 (fi + kS) and satisfies : 

V(t, x) e Ui + «B, a < |V^i(*, x)| < A (Al) 

V(t,x)efii + KB, |^(t,x)|</3, (A2) 

V(t, x) G + k5, |D 2 ft(t,x)|<M (A3) 

and 

V(M)efi» + n5, fe(t,x)| + |a t V^(t,x)| < M. (A4) 
For all t £ I, we denote by 

I{t,x) := {i, gi(t,x) = 0} 

the set of "active contraints" and for p > 

I p {t,x) := {i, gi(t,x) < p} . 

We suppose that there exist constants p, 7 > such that for all x £ Q{t) and all nonnegative 
reals A,- 



\i\V gi (t,x)\ < 7 

ie/ p (t,x) 



XiVgi(t,x) 

i£l p {t,x) 



Theorem 1^.1. — Under the assumptions (Al), (A3) and (Rq), there exists r] := r/(a, M, 7) 
swc/i i/ia£ the set Q(t) is rj-prox-regular for all t £ X. 

Proof: The time t is fixed in this proof so for simplicity we omit it in the notations. We will 
follow the arguments and the ideas of [42] (Subsections 3.1 and 3.2) and [21] (Subsection 2.2) 
where the desired result is already proved in the case of convex constraints gi. 
First let us study the set Qi for a fixed index i. We refer the reader to Proposition 3.2 of [42] 
for the following well-known fact. Due to the assumptions, Qi has a C 1 -boundary 

dQi = {ie R d , 9i (x) = 0} 

and for x £ dQi, its proximal normal cone is given by 

N(Qi,x) = -R+V gi (x). 

We now want to check that Qi is uniformly prox-regular with a constant r/o- It suffices to check 
the hypomonotonicity property: for all x £ dQi and y £ Qi 

(x - y, -Vgi(x)) > - — \x - y\ 2 \Vgi(x)\. (65) 
2% 
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Indeed since gi{x) = 0, a first order expansion together with Assumption (-43) give 

... M l2 

< gi{y) < {x-y, -V^(x)) + —\x - y\ . 

So (65) is satisfied with r]Q := a/M, hence Qi is 770-prox-regular (thanks to Proposition 2.6). 
Then let us study Q the intersection of sets Qi. We first have to prove that for all x G dQ 

N{Q,x)= N(Qi,x) = - R+ V9i(x). (66) 

i£l(x) iel(x) 

Let us denote 

M x •■= - Yl ^ + Vgi(x). 

The inclusion N x C N(Q, x) is proved in Proposition 2.16 [21] (this part did not use the convexity 
of the functions gi) and so we just deal with the other one. By the way, we point out that M x 
is the polar cone of 



T* := \z G 



Vi€/(x), (Vgi(x),z) >o}. 



So using the orthogonal decomposition related to polar cones (see [24]), any v G N(Q,x) can be 
written v = w + z = Pj\f x v + Pt^, with wJ-z. Suppose z 7^ 0. Since v G N(Q,x), there exists 
t > such that x G Pq(x + tu). Let 

9«( x ) 

s = mm(t, r) with r := mm 



i£/(a;) (2/3 + 5a)\z\ ' 

by the well-known property of the projection, the inequality s < t implies 

x G Pq{x + sv). (67) 

From Lemma 5.2 in [6], we know that Q satisfies the second property of the admissibility: there 
exist a bounded covering of dQ with balls B((x p , r)) p , a collection of "good direction (tt p ) p " and 
constants p, 5 such that for all x G B(x p , 2r) and all v G N(Q, x) 

(u,u p ) > <5M- 

There exists p such that x G B(x p , r) and we set 

x := x + sv — sw — es\z\u p = x + sz — es\z\u p , 

where e will be later chosen small enough (e << 1). We claim that x G Q. Indeed thanks to 
Assumption (A3) we have for all i, 

gi(x) > gi(x) + s(Vgi(x),z - e\z\u p ) - s 2 M(l + e) 2 \z\ 2 . 

Consequently, 

Vi G I(x), gi{x) > se\z\5a - s 2 M(l + e) 2 \z\ 2 > 



if s < 



eSa 



\z\M(l+e)*- 

CJ<i (x^) 

Furthermore, if i 4 I(x), then s < r < -n. 7 . , • Hence 

^ K ' ~ ~ (eSa + (1 + e)0)\z\ 

9i(x)>gi(x)-8P\z\(l + e)>0. 
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That is why x G Q as soon as s < \ z \M^i +t ^ an d i n this case cIq(x + si;) < \x + sv — x\ 



w + eztt 



:= syA. By expanding 



.4 



w + e b ki r 



= | it? | +e \z\ +2e\z\(w,u p ) 
< \w\ 2 + e 2 \z\ 2 + 2e|2|M 



< \w\ z + -|ar| 



< \v\ 



I 1 2 i I 1 2 
|w| + |z| , 



for e < i- 



r. Finally we have obtained that for small enough e and s 



2 4H+M 

c?q(x + sv) < s\v\ 

which leads to a contradiction with (67). So we conclude that z = 0, which completes the proof 
of (66). 

Finally, the prox-regularity of the set Q is shown by invoking the "reverse triangle inequality" 
Assumption (R p ), as done in Proposition 2.17 [21]. □ 

Remark Jf.,2. — Note that for all t G X, x G Q(i) 

T t , x := {z G K d , Vt G /(t,a;), (V^(t,x),z) > o] = T Q(t) (x) 
Indeed T^ x and Tqu\(x) are two convex cones whose polar cones are equal: 

i£l(t,x) 

Proposition 4-3. — Let Q be the set-valued map defined at the beginning of the current section. 
The set O := {(t, x) G X x WL d , x G Q(t)} is uniformly prox-regular. Moreover for all (t, x) G 

o 

with t G X, the set Ct jX (defined in Definition 1.1) verifies 

C t , x = [z G R d , Vi G /(i,z), %<(i,x) + {z,V x9i (t,x)) > o} . 

Proof: The set is given by the functions hi as 

n ■= |(t jX ) G R^ 1 , Vi G {0,...,p + l}, &i(t,x) > 0, } 

with /ij := gi for i G {1, --,p}, ho(t,x) := t and h p +i(t,x) := T — t. In order to apply Theorem 
4.1, we check the different assumptions relatively to hi', for all i G {0, ...,p + 1} 

V(£, x) G rij + K.B, min{a, 1} < |V(£ ja .)/ij(i, x)| < max{2/3, 1}, 



V(t,s) G Sli + KB, \D 2 ^ x) hi(t,x)\ < 2M, 



and 



A ! |V M / il (i,x)|< 7 (l + ^j 



iei(t,x) 



iei(t,x) 



(69) 



(70) 



Indeed Properties (68) and (69) are obvious for i = and i = p + 1. Let z G {l,...,p} and 
(t,x) G fij + K-B. Inequality (68) is proved by 

a < \V x9i (t,x)\ < \V {t , x) hi(t,x)\ < \V x9i (t,x)\ + \d t9i {t,x)\ < 2/3. 

The second one (69) is due to Assumptions (A3) and (-A4): 



\D 



l, x) hi(t,x)\ < (\d 2 9i (t,x)\ 2 + \D 2 x9i (t,x)\ 2 + 2\d t V x hi(t, x)\ 2 f /2 < 2M 
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Concerning Assumption (70), we just deal with the case where I(t,x) n {1, 7^ (else the 
inequality is obvious). For all i G {1, ...,p}, we have \dtQi{t,x)\ < ^|V X t/j(i, x)| which implies 

\V it>x) hi(t,x)\ < {\d t gi(t,x)\ 2 + \V xgi (t,x)\ 2 ) 1/2 < (l + ^) \V xgi (t,x)\. 

For i = 0ori=p+l, it comes similarly 

\^(t,x)hi(t,x)\ = 1 < -\V x g io (t,x)\ 

where io G I(t,x) n {1, Then (70) is involved by Assumption (R p ). By Theorem 4.1, we 

conclude to the uniform prox-regularity of Q. 

o 

Thanks to Proposition 2.11, we know that for all (t,x) G £1 with t£l 

Ct,x = {ue R d , G Tg((t,x))} = <ju G M d , (l,u) G T a ((t,x))} 

because O is uniformly prox-regular (see Corollary 6.30 [36]). By Remark 4.2, we get 
T n ((t, x)) = e Vi G J(t, x), (V (t|a .) /ii(i, x),z)> o} 

= {z = (zi,z u ) G Rx R d , Vi G J(i,x)n{l,...,p}, «i%i(t,x) + (V x y i (t J x),z u ) > o} , 

o 

since t G X. Consequently, u G Cj >a; if and only if dtgi(t,x) + {u,V x 9i(t,x)) > for all i G 
I(t,x) = /(t,x)n{l,...,p}. □ 

Proposition 4-4- — Under the previous assumptions, the set-valued map Q(-) is admissible. 
More precisely, for all (io>9o) G ^ £/iere exist a "500c? direction u" and constants r,S such that 
for all (t,q) G QnB((to,qo),2r) and all proximal vector v G N(Q(i), g) = — ^ igJ ^ ^ 1R + V5i(t, (?) 

(u,u)>*H. (71) 
Moreover the set-valued map Q is Lipschitz continuous on I. 



We refer the reader to Lemma 5.2 in [6] for the admissibility property and to Proposition 2.11 
of [43] for the Lipschitz continuity. 

The case of a set-valued map Q has already been studied in [5]. We look for explaining that the 
current results of existence solutions for differential inclusions covers the one obtained in [5] . 

Proposition 4-5. — Under the above assumptions, Problem (9) is equivalent to the following 
one : find q G VF 1,00 (X, R d ) , q G BV(l,R d ) and time-measures Aj G M + (l) such that 

' Vi G X, q(t) G Q(t) 

v 

dq = f(t, q)dt + ^ V g <?)<iA; 
i=l 

< s<w(Ai) C {* , gi(t, q(t)) = 0} for all i (72) 
Vtel, q(t+)=P CtMt) q(t-) 

9(0) = 90 
> 9(0) = u . 

Proof: By the characterization of the proximal normal cones (see Remark 4.2), it is obvious 
that a solution of Problem (72) is a solution of Problem (9) too. Indeed, consider q, Ai,...,A p 
a solution of Problem (72). Then we define the measure A = Ai + ... + A p in order that each 
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measure Aj is absolutely continuous with respect to A. There also exist bounded nonnegative 
and mesurable functions l{ such that for all i = 1, ...,p: dXi = lidX. Then 

dq = f(t,q)dt + ^2V q g i (t,q)dX l = f(t,q)dt+ ( ^ V, 9i(t, q)h{t) J dA. 

i=l \t=l / 

By writing £(i) = >Xi V 9 ft (t, q)h (t) G -N(Q(t),g(t)), we set £(t) = * &*) is 

non- vanishing and £(t) = else. So, we obtain 

d<7 = f(t,q)dt + €(t)dv(t) 

with d^ = |£(£)|dA. Since the functions £i are bounded by 1, the function £ is bounded on I and 
f G .M+(X). We conclude that q is a solution of Problem (9). 

Let us explain the other relation. Let q be a solution of (9). By definition, there exist a function 
k G £>y(Z, R d ) and a measurable map £ : X — > M. d such that 

dq + dk = f(t, q)dt 

and ^ 

k(t)= [ t(s)d\k\(s) 







with £(s) G N(C(s), and |£(s)| = 1. The technical difficulty is to represent the map £ in 
terms of Kuhn- Tucker multipliers (this corresponds to a problem of selection for a particular 
set-valued map). For all t G X, we define 



F{t) : = j(Ai, Ap) G £(*) = - £ A;V ? ft(t, 9(0), A; / only if 9i (t, q(t)) = J . 

We are now looking to obtain a measurable selection of the set-valued map J- . It is easy to see 
that T takes non-empty closed values. Moreover its graph Tjr 

ry : = {(t,Ai,...,A p )GXx (R+y, (a 1; ...,a p ) e?(t)} 

is given by 

:=g-\{o})n{ix(R+)P) 

where Q is defined as follows 

g(t,X 1 ,...,X p ) :=£(f) + ^ XiV q gi(t,q(t)) + ^ |Ai|xj(t,q(t))=(i) 

iel(t,q(t)) \i0(t,q(t)) 

where e is any vector e G M d \ {0} and Xi(t,q(t)) c is the- characteristic function: 

m _/0, i€l(t,q(t)) 
Xl(t, q (t)<V ~ | +OC) i£i(t iq {t)). 

Consequently, since £ is measurable and t — > I(t, q(t)) is upper semicontinuous, we deduce that Q 
is measurable and so it follows that the graph Tjr is a measurable set. Then, Theorem 8.1.4 in [1] 
yields that the set-valued map J- is measurable and that J- admits a measurable selection. Let 
us write (oi, a p ) for such a measurable selection of T . Then, we let the reader to check that 
the measures dXi := etjd|A;| are solutions of Problem (72) (since the reverse triangle inequality 
(R p )i it comes that the functions a, belong to L°°(I,M)). □ 

We have checked that the set-valued map Q takes uniformly prox-regular values (see Theorem 
4.1). Moreover the Lipschitz regularity and the admissibility property have already been proved 
(see Proposition 4.4). We can now apply Theorem 1.3 and we get the following one: 

Theorem 1^.6. — Under the above assumptions ((Al)-(R p )) and with f satisfying (12) and 
(13), the problems (5) and (72) are equivalent and admit at least one solution on any time 
interval X. 
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We would like to finish this section with a local existence result for a constant set Q(t) 
defined by time- independent constraints. 



Q 



Theorem 4-7. — Let k > and consider Q := Hf=i { x £ 9i( x ) — 0}> where gi G C 2 (Q + 
kB(0, 1)) satisfy (Al ), (A3) and (Rq). Then for all initial data go E Int(Q) and uq E i/iere 
exisi To := To(|i*o|) > and q solution of (2) on [0, To] where f satisfies (12) and (13) with 
F E L°°(X). 

Proof: First, Q is uniformly prox-regular by Theorem 4.1. Thanks to Theorem 5.1, there 
exists To := To(|«o|) > such that the computed velocities (u^jh & re uniformly bounded in 
L°°([0, To], M. d ) (for h going to 0). We can also conclude that all the computed solutions qh are 
uniformly bounded: there exists L > such that for all small enough h > and all t E [0, To] : 

Due to Lemma 4.8, there exist p := p(L) and 7 := 7(T) > such that (R p ) holds for every 
q E B(0, 2L). By Lemma 5.2 in [6], the set Q verifies a "local admissibility" property: there exist 
5, r > 0, sequences (x p ) p and with \u p \ = 1 and x p £ Q H B(0,2L) such that (B(x p ,r)) p is 

a bounded covering of the boundary dQ n i?(0, 2L) and 

Vp, Vx E dQ n i?(0, 2L) n B(x p , 2r), WeN(Q,x), (v, u p ) > S\v\. 

We let the reader to check that the proof of Theorem 1.3 still holds since we consider points 
q% E Q n 5(0, L) (indeed we can choose r < L in order that for all g? E -B(0, T), the whole ball 
B(qf t ,r) is included in B(0,2L) where the admissibility property is verified). Thus we get the 
existence of a solution on [0,To]. □ 

Lemma 4.8. — Let L > then there exist p := p(L) and 7 := 7(T) > such that (R p ) holds 
for every q E B(0, 2L). 

Proof: Assume that (R p ) does not hold for all p > in B(0,2L). So there exist a sequence 
(<7n)n E B(0,2L) and (Xf)i, n such that 



n 



E A ? v ^«: 



E A riv ffi ((?„)i > 

ie/i/„(gn) 

with A" > 0. Without loss of generality, we can assume that Sie/ 1/n (g„) 

By compactness, up to a subsequence, there exist g E i?(0, 2L) and (Ai)i=i ) ... )P such that q n 
and for all i E {1, ...,£>}, A" — > Aj. It comes with (A3) 



n 



So 



E A " v ^) 

*e^l/n(?n) 



E A " v ^) 

«eJ"i/„(g n ) 



n 



M|g n - g| < E A "l V ^(<z)l + M|g n - q\ 



<- E A ?|Vffi(9)|+fM + -)|g„-q 



< - + M+— g re -g. 
n \ n / 

Furthermore, since E C 1 (M rf ), gj is uniformly continuous and so it is easy to check that 
Li/ n (q n ) C /(g) for n large enough. For such n, (Rq) and (Al) imply 



E A nv^)i< 



E A " v ^(<?) 

iell/„(?n) 
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Finally, we obtain 



a 



< - + I 



7 re 

which leads to a contradiction for re — > oo. 



re J 



□ 



In Theorem 4.7, the assumptions on the gradients (Vpj)j are weakened with respect to the 
existing results (see the introduction). Indeed the gradients of active constraints are usually 
supposed to be linearly independent, which implies that Assumption {Rq) holds on any com- 
pact (see Remark 4.9). Assumption {Rq) permits also to deal with a large number of active 
constraints, while the linearly independence requires that \I{q)\ < d. In addition, we emphasize 
that the classical geometrical assumption: for all i ^ i' G I{q) 

(Vg t {q),Vgii{q)) <0 

is not required. 

Remark 4-9. — Let K Cl d be a compact set and assume that for all q G QCiK, (V ' gi{q))i^i( q ) 
is linearly independent. This collection can be extended to a basis ofM. d , denoted by (e?)j = i ... d- 



Let us consider 



d the associated dual basis such that for all i,j G {1, .., d} 



4) 



1, i = j 
0, i^j. 



Then for q G Q n K and nonnegative reals {^i)iel(q)> & comes 

d 

E *iiv*(fl)i <p E ^f(Vft(?))</3 E EM^W) 



*6/(g) 



*el(g) i=l 



ie/(g) 

d 

^^E £ l ( E A ^(?) ] <p 



J'=l 



E AiVpi(g) 
*ei(?) 



5o 3 we obtain (Rq) on K with 



7 := 



/3 sup 



3=1 



< OO 



fry Lemma 7 of [32] . 



5. Appendix 

This section is devoted to another proof of uniform bounds for the computed velocities in 
L°°{l,R d ) (Proposition 3.2). However, we are looking for a proof in the framework of a constant 
set C, requiring only the prox-regularity property on C and we do not assume it is admissible. 
We will see that we can just obtain a local result: the velocities are bounded on a time interval 
[0, To], where the time To depends on the initial conditions. It is not clear how can we extend 
the proof to a global result without extra properties such the admissibility. 

Theorem 5.1. — Assume that C(-) := C where C is a uniformly prox-regular (possibly not 
admissible) set and f satisfies (12) and (13) with F G L°°(I). Then, the computed velocities 
(still by using the scheme (26)) {uh)h<i o,re uniformly bounded in L°°{I,M. d ) on [0, To], where 
To depends on the initial data uq as follows: 



2{J + 1) (2\u \ + 3HFIU + ^Fj^) ' 
with a numerical constant J (defined in Lemma 5.3). 
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We first detail three technical lemmas and postpone the proof after them. 
Lemma 5.2. — Let n be an integer and 7 e]0, 1/4] (later defined in Lemma 5.3). Then if 
K +2 -Q% +1 \<V and |g™ +1 - q n h + h 2 f% +1 \ < 7 r? 

then 

K +2 - q n h +1 \ < K +l -qt + h 2 fZ +1 \ + J\q n h +l - q n h + h 2 f% +1 \ 2 . 

Proof: By definition, for each integer n, 

q n h +2 := P C [2q n h +1 - q n h + tfff 1 ] . 

It follows that with z := 2q^ +1 — q% + h 2 f£ +1 , z — q^ +2 is a proximal normal vector of C at the 
point q^ +2 - The hypomonotonicity property of the proximal normal cone yields 



It also comes 

n+2 „n+l|2 



n+l „n+2 „n+l 



1 



n+2 „n+l|2 



- n+1 ' + \q n h +2 -q n h +1 \). 



2t] 



Using Cauchy-Schwarz inequality, we deduce 



K +2 - € 



1 



\*-<fir L \<2j\<c-ti 



Thanks to Lemma 5.3 (later proved) with a :- 
comes 

K +2 -€ +1 \<\*-€ +1 \ + J\*-€ +1 \ 2 



z-q n h +l \ + K +2 -q n h +1 \). 
n+2 _ ? n+l| < v and b:= \ z _ < ^ it 



□ 



Lemma 5.3. — Let a and b two nonnegative reals satisfying 

a 2 + (b- 2r])a + 2i]b > 0. (73) 
Then for some numerical constant 7 G]0, l/4[ and J > (independent on a, b), we have 



a < rj 
b < 777 



a < b + Jb 2 . 



(74) 



Proof: We remark that (73) is a second degree polynomial function with respect to a, whose 
discriminant is given by 

A := (b - 2r]) 2 - 8i]b. 

The real A is nonnegative as soon as b < 777 < 77/4. Then we know that (73) implies 

a x + [ 

with 



It comes 



-6 + 277 VA 
x ± := 7, ± — — . 



A = 4t7 2 



b b 2 

1 " 3 - + TT 
77 4?7 Z 



Hence, 



< M6 2 



for some numerical constant M = M( 7 , 77). 

6 



+ 77 ± 77 



36 
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So X- < b+Mb 2 /2 and > 2r/— 2b— Mb 2 /2. Furthermore, there exists a small enough constant 
7 > such that 2777 + M^rf) 2 /2 < 77, which implies with b < 777 that > r/. Consequently, 
since a and a < 77, we deduce that necessarily 

M. 



a < x_ < b + ^6 2 , 



which concludes the proof of (74). □ 
Lemma 5.4- — Let J' be a fixed positive real. We denote by (j) n j, the n-th iterated of 

4>ji \ = x 1 y x ~\~ J'x 2 . 
For every x > with J'nx 7^ 1, we have 

<fi,(x) < x ( t J ' n f ~1=± ) = x h+ {J'nx) + ■■■ + (J'nx) 2 "- 1 ) . (75) 
\ J'nx — 1 J ' 

Proof: It is obvious that (jfj, is a polynomial function of degree 2 n with nonnegative coefficients. 
So we know that (jfj, can be written as 

fc=i 

We want to prove that for all n > 1 and A; £ {1, 2 n } 

4 n) < (J'nf- 1 . (76) 

It is obvious that (76) holds for n = 1 (in fact, there is equality). Let us assume that (76) holds 
for an integer n and prove it for n + 1. Since 

r/ l (x) = rAx) + J'rAx) 2 , 

for all k G {l,...,2 n+1 } 

4 n+1) = 4 n) + ^E«r4- ) J 

fe-i 
i=i 

< (J'n) fe_1 + J'(k-l)(J'n) k - 2 

< {J'{n + l)) k -\ 

(n) 

By agreement, a k is set equal to if k > 2 n . 

That ends the recursive proof. □ 

We now come back to the proof of Theorem 5.1. 

Proof of Theorem 5.1: We define a sequence (x n ) n as follows: 

x° := \ql-ql\+h 2 \\F\\^ + h^WU and x n+1 = 4>J + i(x n ), 

where J is introduced in Lemma 5.3 and functions <pj in Lemma 5.4. We choose h < min{To, 1} 
small enough in order that xo < ^7?? (this is possible since \q\ — q^\ < h\u l h \ < 2K\uq\ + 2/i 2 ||F|| 0O , 
see Lemma 3.1). 
Now we set 

V := {n < T/h, \q n h +l -q n h \<n and \q n h - g^l + h 2 |/£| < 7 n} 

and 

P := min{n > 0, n <£ V} - 1. 
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We want to prove that P = N where N := Tq/H. Let fix n £ {1, P}, since \q^ +l — qV\ < rj 
and IgJJ — + ^- 2 |/^| < 7Vi thanks to Lemma 5.2 it comes 

K +1 - tf| < |tf - qt 1 + h 2 fjl\ + J\ql- q^ 1 + h 2 f%\ 2 
<d>.,{\q n h -q n h - l + h 2 fZ\). 

Consequently 

Wh +1 - qt\ + h 2 \ft +l \ < M\q n h - q n h ~ l \ + h 2 \m) + tfllFlU 

<<l>j{\q n h -q n h - l \+h 2 \m) + (x'f 

<M\q n k-q n h - 1 \ + h 2 \m) + ^ n ~ 1 ) 2 , 

where we have used that the sequence (x n ) n is increasing. Consequently, we easily deduce by 
iteration that for all n = 1, P, 

\q n h +l -q n h \+h 2 \fZ +l \<x\ (77) 
Indeed, (77) is satisfied for n = 0. Moreover assuming (77) for n — 1, we have 

k +1 - Q n h \ + h 2 \K +1 \ < <mk - ^r 1 ! + ^u n i) + (^ n - 1 ) 2 < ^j+i^- 1 ) = 

because <f>j is nondecreasing on [0,oo[. 

Thanks to Lemma 5.4 (with J' = J + 1), we have for every n < P 
Moreover for all n < N = To/h 



(J + l)nx° < (J + 1)^ (|<£ - gg| + h 2 ||F|U + 



<(J+1)T (|«i| + /i||F|| 00 + #1 



< (J + 1)T [2\uo\ + 3/iHFHoo + VFT 
1 

where we have used Lemma 3.1 to estimate \ul\ by 2|uo| + 2/illFjloo and the definition of To. 
Thus for n = P, we get 

k P+1 -^l+/ l 2 |4 P+1 |<2x°< 7 r ? . 

Due to Lemma 3.1, it comes 

\lh +2 ~ < 2|< +1 - 9 P | + 2h 2 \f^\ < 2 7r] < r,. 

Finally, P + 1 £ "P and so P = N = To/h. We also conclude that 

sup |u]J | < sup ~ 9/1 



0<t™<T o n<7V |^| 

< 2— 
h 

<2f|4| + /t||F|| 00 + v ^| 



< 2 ^2|u | + 3/i||F||oo + vIM, 
which is uniformly bounded with respect to h when h goes to 0. □ 
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